Energy  Budget  of  Nonlinear  Internal  Waves  near  Dongsha 


Ren-Chieh  Lien 
Applied  Physies  Laboratory 
University  of  Washington 
Seattle,  Washington  98105 

phone;  (206)  685-1079  fax;  (206)  543-6785  email;  lien@apl.washington.edu 

Erie  A.  D’Asaro 
Applied  Physies  Laboratory 
University  of  Washington 
Seattle,  Washington  98105 

phone;  (206)  685-2982  fax;  (206)  543-6785  email;  dasaro@apl.washington.edu 
Award  Number;  N00014-05-1-0284 


LONG-TERM  GOALS 

Our  long-term  seientifie  goal  is  to  understand  the  meehanisms  by  which  mixing  occurs  in  the  ocean  and 
thereby  help  develop  improved  parameterizations  of  mixing  for  ocean  models.  Mixing  within  the 
stratified  ocean  is  our  particular  focus  as  the  complex  interplay  of  internal  waves  from  a  variety  of 
sources  and  turbulence  makes  this  a  current  locus  of  uncertainty.  In  this  study,  our  broad  focus  is  on  the 
energy  sources  of  nonlinear  internal  waves  (NLIWs)  in  a  complex  environment  of  strong  internal  tides 
and  abrupt  topography  (continental  shelf  and  slope).  We  expect  a  rapid  evolution  of  internal  tides  and 
NLIWs,  and  we  aim  to  understand  their  dynamics,  energy  cascade,  and  role  in  mixing. 

OBJECTIVES 

The  primary  objectives  of  the  present  project  are  1)  to  identify  the  generation  site  and  understand  the 
generation  mechanism  of  NLIWs,  2)  to  understand  the  evolution  of  NLIWs  interacting  with  abrupt 
topography,  3)  to  quantify  the  energy  budget  and  energy  cascade  from  internal  tides  to  NLIWs,  and  4)  to 
quantify  the  seasonal  variation  of  the  energy  of  NLIWs  near  Dongsha  Island  in  the  northern  South  China 
Sea  (SCS).  Our  particular  interest  is  to  understand  the  energy  cascade  from  barotropic  tides,  internal 
tides,  NLIWs,  to  turbulence  mixing  in  the  northern  SCS,  and  to  understand  the  evolution  of  NLIWs 
interacting  with  the  shoaling  continental  slope. 

APPROACH 

Our  approach  is  to  take  direct  observations  of  NLIWs  near  Dongsha  Island  where  most  of  NLIW 
activities  are  captured  in  satellite  images.  In  April-May  2005,  we  conducted  a  two-week  observational 
experiment  near  Dongsha.  Large-amplitude  NLIWs,  greater  than  150  m,  and  strong  turbulence  mixing 
were  observed  by  the  Lagrangian  float  and  shipboard  sensors  including  ADCP,  CTD,  EK500,  and 
X-band  marine  radar.  Combined  remote  sensing  and  in-situ  measurements  provided  detailed  properties 
of  large-amplitude  NLIWs.  In  June  2006,  an  array  of  ADCP  moorings  was  deployed  along  the  prevailing 
path  of  NLIWs  near  Dongsha.  Long-term  observations  of  NLIWs  will  be  obtained  allowing  us  to  (1) 
quantify  the  seasonal  variation  of  NLIW  energy,  (2)  map  the  geographical  distribution,  (3)  better 
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understand  the  dynamies  of  NLIW  evolution  over  the  shoaling  topography,  and  (4)  assess  the  model 
predietion  skill  of  NLIWs. 

WORK  COMPLETED 

A  set  of  mooring  ADCP  data  provided  by  Dr.  David  Tang  of  National  Taiwan  University  was  analyzed 
and  results  were  published  (Lien  et  al,  2005  and  Chang  et  ah,  2006).  These  ADCP  measurements 
showed  strong  and  highly  predietable  NLIWs  on  the  eontinental  slope  and  plateau  in  the  vieinity  of 
Dongsha  Island.  Strong  energy  flux  divergenee  of  NLIWs  and  internal  tides  are  found  aeross  the 
Dongsha  plateau. 

Measurements  taken  in  April  18-May  1  of  2005  are  analyzed.  Dr.  Sue  Moore,  the  whale  expert  at  APL, 
was  attraeted  to  our  video  eapture  of  a  sehool  of  pilot  whales  in  the  NLIWs.  A  paper  is  to  be  published  in 
Marine  Mammal  Seienee  (Moore  and  Lien,  2006).  We  foeus  our  analysis  on  deseribing  1)  the  rapid 
evolution  of  NLIWs  on  the  shoaling  eontinental  slope,  and  2)  the  surfaee  signature  indueed  by  NLIWs. 
Results  of  these  analyses  were  presented  in  the  2006  Oeean  Seienees  Meeting  held  in  Honolulu  and  will 
be  deseribed  in  the  next  seetion. 

A  set  of  three  bottom-mounted  Long  Ranger  ADCPs  were  deployed  in  June  2006.  These  moorings  are 
seheduled  to  take  more  than  one  year  of  measurements  of  NLIWs.  They  will  be  servieed  every  3-4 
months  for  retrieving  data,  replaeing  batteries,  and  maintenanee. 

RESULTS 

Energy  Flux  and  Flux  Divergence  across  the  Dongsha  Plateau 

Three  sets  of  ADCP  measurements  taken  on  the  Dongsha  plateau,  on  the  shallow  eontinental  shelf,  and 
on  the  steep  eontinental  slope  in  the  northern  SCS  are  analyzed  (Chang  et  ah,  2006).  The  data  show 
strong  divergenees  of  energy  and  energy  flux  of  NLIWs  along  and  aeross  waves’  prevailing  westward 
propagation  path  (Fig.  1).  The  NLIW  energy  flux  is  8.5  kW  m-1  on  the  plateau,  only  0.25  kW  m-1  on  the 
eontinental  shelf  220  km  westward  along  the  propagation  path,  and  only  1  kW  m-1  on  the  eontinental 
slope  120  km  northward  aeross  the  propagation  path.  Along  the  wave  path  on  the  plateau,  the  average 
energy  flux  divergenee  of  NLIWs  is  ~0.04  W  m-2,  whieh  eorresponds  to  a  dissipation  rate  of  0(10-7 
-10-6)  W  kg-1  for  NLIWs.  Combining  these  observations  and  model  results,  we  eonelude  a  seenario  of 
NLIW  energy  flux  in  the  SCS.  NLIWs  are  generated  east  of  the  plateau,  propagate  predominantly 
westward  aeross  the  plateau  along  a  beam  of  -'100  km  width  that  is  eentered  at  -'21oN,  and  dissipate 
nearly  all  their  energy  before  reaehing  the  eontinental  shelf. 

Properties  of  NLIWs  on  Shoaling  Continental  Slope 

In  April-May  2005,  our  measurements  were  taken  mostly  along  2  UN  between  1 18  on  the  upper 
eontinental  slope  and  1 16.5°E  on  Dongsha  plateau  of  about  350-500-m  depth  (Fig.  2)  using  shipboard 
ADCP,  CTD,  EK500,  marine  Radar,  XBT,  and  a  Eagrangian  float.  The  float  measures  internal  waves 
and  turbulenee  mixing. 
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The  NLIW  evolves  rapidly  as  it  propagates  westward  up  the  eontinental  slope  onto  Dongsha  plateau.  On 
April  28,  the  details  of  the  evolution  were  eaptured  by  shipboard  observations  and  summarized  in  Fig.  3. 
The  NLIW  first  appears  as  a  single-depression  wave  with  amplitude  >  170m.  Within  two  hours  and  in  a 
distanee  of  10  km,  a  small  wave  is  born  trailing  immediately  behind  the  original  wave,  identified  by 
marine  radar  measurements.  The  small  wave  propagates  slower  than  the  leading  large  wave,  and  is 
further  separated  from  the  leading  wave  as  they  propagate  westward  aeross  Dongsha  plateau. 

The  dynamies  of  this  fission  proeess  is  important  for  the  energy  evolution  of  NLIWs,  but  is  yet  not  well 
known.  A  eloser  look  of  the  veloeity  strueture  and  the  speed  of  NLIWs  provides  a  hint  of  the  evolution 
(Fig.  4).  As  the  NLIW  shoals,  the  wave  speed  deereases  immediately  responding  to  the  loeal  depth, 
whereas  the  maximum  fluid  veloeity  remains  unehanged  for  a  while.  A  trapped  eore  is  formed.  Within 
the  trapped  eore,  Ri  is  smaller  than  ki,  and  the  dissipation  rate  of  turbulenee  kinetie  energy  is  O  (10-5)  W 
kg-1  observed  by  the  Lagrangian  float  (Fig.  5).  An  overturning  seale  of  75  m  is  observed.  More  than  3.5 
hours  after  the  NLIW  reduees  its  speed,  the  maximum  horizontal  veloeity  of  the  NLIW  begins  to 
deerease  gradually,  presumably  as  a  result  of  the  eombined  fission  proeess  and  dissipation.  About 
another  3  hrs  later,  the  trapped  eore  disappears  (Fig.  4).  The  eartoon  in  Fig  6  illustrates  the  proposed 
seenario  of  the  trapped  eore  formation,  fission,  dissipation,  and  evolution  of  NLIW  on  the  shoaling 
topography.  The  hypothesis  needs  to  be  further  investigated  by  observations  and  model  simulations. 

Remote  Sensing  and  Prediction  of  Properties  of  Nonlinear  Internal  Waves 

We  analyze  the  simultaneous  set  of  marine  radar,  shipboard  ADCP,  CTD,  and  eeho  sounder 
observations  of  large-amplitude  NLIWs  in  the  SCS.  Primary  eauses  of  strong  surfaee  seattering  in  radar 
observations  are  winds  and  NLIWs.  In  the  absenee  of  NLIWs,  the  wind  speed  is  strongly  eorrelated  with 
the  radar  seattering  intensity  with  a  eorrelation  eoeffieient  of  0.87,  but  the  eorrelation  disappears  when 
NLIWs  are  present.  We  analyze  properties  of  NLIWs  and  radar  surfaee  seattering  intensity  during  events 
of  strong  NLIWs  and  eonstruet  a  eomposite  strueture  of  NLIWs  and  the  eorresponding  surfaee  seattering 
strength. 


The  average  half  width  (  of  NLIWs,  i.e.,  the  spatial  seale  of  one  half  of  the  maximum  amplitude  of 

NLIWs,  is  1.09  ±  0.2  km  (Fig.  7).  The  average  half  width  of  the  seattering  intensity  and  strain  rate  is 

0.57  ±  0.12  (top  panels  in  Fig.  7)  and  0.46  ±  0.06  (middle  panels  in  Fig.  7),  respeetively.  The 

A 

strong  surfaee  seattering  leads  the  maximum  amplitude  of  NLIWs  by  0.46  ±  O.I  I  .  Surfaee 

seattering  intensity  measured  by  the  radar  looking  from  behind  the  NLIW  (blue  eurves  in  the  top  panels) 
is  76%  of  that  measured  by  the  radar  looking  from  ahead  the  NLIW  (red  eurves  in  the  top  panels).  In  the 
du 

divergenee  zone  (  dx  >o),  the  rear  portion  of  NLIWs,  surfaee  seattering  intensity  observed  by  the  radar  is 

du 

not  signifieantly  different  from  the  baekground  (Fig.  8).  In  the  eonvergenee  zone  (  dx  <o),  the  front 
portion  of  NLIWs,  the  surfaee  seattering  intensity  is  positively  proportional  to  the  magnitude  of  the 

du 

strain  rate.  A  linear  regression  line  fits  reasonably  well  in  -1.8x10-3  s-1  <  dx  §.1  ^vvtth  a  eorrelation 
eoeffieient  -0.93  (red  eurves).  Over  the  entire  strain  rate  observed,  an  aretangent  model  fits  better  (blue 
eurves).  Our  results  eonelude  that  strain  rates,  horizontal  seales,  and  positions  of  NLIWs  are  predietable 
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by  remote  sensing  observations.  The  applieation  of  these  results  on  satellite  derived  surfaee  signatures  of 
NLIWs  is  praotieal. 


IMPACT/APPLICATION 

Our  analysis  eoneludes  that  NLIWs  evolve  rapidly  aeross  the  upper  flank  of  the  eontinental  slope  and  the 
Dongsha  plateau  via  eomplieate  proeesses,  e.g.,  the  formation  of  trapped  eores,  and  the  development  of 
wave  train.  These  proeesses  are  responsible  for  the  strong  dissipation  of  NLIWs  in  SCS.  Our  analysis  of 
eombined  remote  sensing  and  in-situ  measurements  yields  a  model  for  predieting  NLIW  properties 
applieable  to  satellite  observations.  Further  long-term  observations  of  NLIWs  in  the  vieinity  of  Dongsha 
plateau  are  in  progress  and  will  provide  a  better  predietion  of  NLIWs  in  SCS. 
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Figure  1:  Flux  divergence  of  NLIWs  on  Dongsha  plateau.  The  average  dissipation 

AF 

s  = - 

rate  is  calculated  as  ,  where  is  the  difference  between  the  vertically 

integrated  energy  flux  of  NLIWs  on  the  eastern  edge  of  the  Dongsha  plateau,  8.5  kW 
m-1,  and  that  on  the  continental  shelf,  0.25  kW m-1,  and  ^  is  the  distance  between 
the  two  stations.  The  flux  divergence  is  ~0.04  W  m-2,  corresponding  to  a  dissipation 
rate  of  0(10-7 )  Wkg-1  for  NLIWs.  The  dashed-line  and  solid-line  boxes  illustrate 
the  horizontal  spreading  of  NLIWs  energy. 
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Figure  2:  (a)  Map  of  the  northern  South  China  Sea,  and  (b)  bottom  bathymetry  along  21°4’N.  In 
panel  (a),  green  curves  represent  surface  signatures  of  NLIWs  identified  in  satellite  images  (Zhao  et 
al.,  2004),  dashed  for  single-depression  waves  and  solid  for  multiple  wave  packets.  The  blue  box 
delineates  the  area  where  multiple  wave  packets  are  mostly  found.  Four  primary  stations  in  our 
April-May  2005  cruise  are  labeled  as  O,  E,  D  and  W.  Shipboard  and  float  measurements  are  taken 
along  0-E-D-W.  Three  magenta  curves  illustrate  isobaric  orientations  on  the  continental  slope.  In 
panel  (b),  two  submarine  ridges  in  the  Luzon  Strait  are  labeled.  They  are  responsible  for  generating 
strong  internal  tides.  Depths  at  four  primary  stations  are  also  labeled.  The  inset  (c)  shows  the 
Lagrangian  float  and  sensors  equipped  on  the  float. 
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Figure  3:  The  top  panel  shows  the  evolution  ofNLIWs  propagating  westward  on  the 
shoaling  continental  slope  captured  by  shipboard  measurements  in  April  28,  2005. 
Solid  dots  represent  the  leading  wave,  and  open  circles  are  trailing  wave.  Amplitudes 
of  NLIWs  are  labeled.  The  bottom  panel  shows  the  bottom  topography.  Two  insets 
are  illustrations  of  radar  images  taken  at  1 745  LT  on  the  upper  flank  of  the 
continental  slope  and  at  2303LT  on  the  plateau.  Red  lines  on  radar  images  indicate 

crests  of  NLIWs. 
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Figure  4:  The  top  panel  summarizes  the  variations  of  the  maximum  horizontal 
velocity  (red  curve)  and  the  wave  speed  (blue  curve)  of  a  large-amplitude  NLIW 
propagating  westward  on  the  shoaling  continental  slope.  The  bottom  panel  shows  the 
bottom  bathymetry.  The  red  box  in  the  inset  delineates 
the  location  of  our  observations. 
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Figure  5:  The  density  contour  obtained  by  the  shipboard  CTD  illustrates  ~75-m 
overturn.  The  Richardson  number  and  e following  the  Lagrangian  float  are 
estimated.  The  red  enclosure  indicates  the  possible  trapped  core. 
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Figure  6:  A  summary  cartoon  illustrates  the  evolution  of  NLIWs  on  the  shoaling 
continental  slope  across  the  Dongsha  plateau. 


Figure  7:  Spatial  structures  of  the  surface  scattering  intensity  observed  by  the 
marine  radar  (top  panels),  the  horizontal  strain  rate  observed  by  the  shipboard 
ADCP  at  14-m  depth  (middle  panel),  and  isopycnal  displacements  at  equilibrium 
depths  of  100  m  and  150  m  (bottom  panels).  The  NLIW propagates  toward  the  left  of 
the  page.  These  properties  are  averaged  over  all  large-amplitude  NLIWI  events.  The 
left  panels  show  the  actual  spatial  scale,  where  X  is  the  horizontal  scale  in  the 
propagating  direction  of  NLIWs.  In  the  right  panels,  X  is  normalized  by  the 
half-amplitude  wave  width  of  NLIWs.  The  surface  scattering  intensity  depends  on 
the  position  of  the  radar  relative  to  the  propagating  direction  of  NLIWs.  The 
intensity  is  stronger  when  the  radar  is  ahead  of  the  NLIW  (red  curves),  and  is  weaker 
when  the  radar  is  behind  the  NLIW  (blue  curves).  Gray  shadings  represent  95% 

confidence  intervals. 
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Figure  8:  Scattering  plots  and  model  fits  between  surface  scattering  intensity 
observed  by  radar  and  the  strain  rate  computed  from  shipboard  ADCP 
measurements.  The  left  panel  shows  observations  (gray  dots)  when  the  radar  is 
ahead  of  the  propagating  NLIWs.  Red  dots  represent  averages  over  constant  grid 
interval,  0.005  s-1,  of  the  strain  rate.  Vertical  lines  represent  95%  error  bars.  Red 
and  blue  curves  represent  the  linear  and  arctangent  fits  to  observations  (red  dots). 

Three  regions  labeled  as  1,  2,  and  3  represent  the  divergence  zone,  the  weak 
convergence  zone,  and  the  strong  convergence  zone.  Symbols,  curves,  and  labels  in 
the  right  panel  are  the  same  as  in  the  left  panel,  but  radar  observations  are  taken 

behind  the  propagating  NLIWs. 
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